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Explicit expressions for the heat and momentum fluxes are given for a low-density multicompo- 
nent mixture in a steady state with temperature and velocity gradients. The results are obtained 
from a formally exact solution of the Gross- Krook model [Phys. Rev. 102, 593 (1956)] of the 
Boltzmann equation for a multicomponent mixture. The transport coefficients (shear viscosity, vis- 
cometric functions, thermal conductivity and a cross coefficient measuring the heat flux orthogonal 
to the thermal gradient) are nonlinear functions of the velocity and temperature gradients and the 
parameters of the mixture (particle masses, concentrations, and force constants). The description 
applies for conditions arbitrarily far from equilibrium and is not restricted to any range of mass 
ratios, molar fractions and/or size ratios. The results show that, in general, the presence of the 
shear flow produces an inhibition in the transport of momentum and energy with respect to that 
O ' of the Navier-Stokes regime. In the particular case of particles mechanically equivalent and in the 

tracer limit, previous results are recovered. 
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^3 I I. INTRODUCTION 

a . 

Needless to say, the study of transport phenomena in fluid mixtures is much more complicated than that of a single 
fluid. Not only is the number of transport coefficients much higher but also they depend on parameters such as the 
mass ratios, the molar fractions and/or the size ratios. These difficulties increase considerably when one attempts 
to analyze far from equilibrium states for which the linear relationships between fluxes and gradients do not apply. 
Due to the the complexity of the general problem, tractable specific situations must be considered. The study of such 
situations (for which a complete description can be given) is of great value since they allow us to gain some insight 
into the understanding of more general problems. 
£ — . In order to capture the relevant aspects of nonlinear transport phenomena, a low-density iV-component mixture 
with short-range interactions can be chosen as a prototype system. In this case, the essential information on the 
mixture is given by the one-particle velocity distribution functions /j(r, v;t) (i = 1, . . . ,N), which obey the set of 
coupled nonlinear Boltzmann equations For states near equilibrium, the constitutive Navier-Stokes equations for 
the heat and momentum fluxes (Fourier's heat law and Newton's viscosity law) can be derived, for instance, by using 
the Chapman- Enskog method H. Nevertheless, beyond the Navier-Stokes domain, the Chapman-Enskog calculations 
are prohibitively difficult and not applicable for large gradients. An interesting problem in which some insight into the 
. mechanisms of heat and momentum transport far from equilibrium can be gained is that of steady planar Couette flow. 
The physical situation of the system is that of a multicomponent gas enclosed between two parallel plates in relative 
motion and maintained, in general, at different temperatures. These boundary conditions lead to temperature and 
O ■ density gradients coexisting with a velocity field. In the steady Couette flow there are two local parameters measuring 
the departure from equilibrium: the shear rate and the thermal gradient. Our aim is to get the relevant transport 
" *~J i coefficients in terms of both gradients, as well as in terms of the parameters characterizing the mixture. 

Unfortunatelly, due to the mathematical difficulties embodied in the nonlinear Couette flow problem, no analytic 
solution to the Boltzmann equation valid for arbitrary values of the shear rate and the thermal gradient is known, even 
for a single gas. Only in the limit of small shear rates, a perturbation solution in powers of the shear field through 
super-Burnett order has been recently worked out for a single gas of Maxwell molecules Q . If one wants to obtain 
the full nonlinear dependence of the transport coefficients on the imposed gradients, either one performs computer 
simulations or on the analytical side one considers kinetic models or the Grad method m. Here, we take the second 
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route and consider a kinetic model. Specifically, we use the nonlinear kinetic model for mixtures proposed by Gross 
and Krook (GK) §. This model is constructed in the same spirit as the well-known Bhatnagar-Gross-Krook (BGK) 
model of a single gas 0, which admits an exact solution for the steady nonlinear Couette flow @J^]. The GK model 
preserves the most important physical properties of the true Boltzmann equation while allowing for more complete 
calculations. 

As mentioned above, the goal of this paper is to obtain exactly the pressure tensor (momentum transport) and the 
heat flux (energy transport) of a multicomponcnt mixture under Couette flow in the context of the GK model. We 
assume that the molar fractions are constant, so that no mutual diffusion appears in the system. On the other hand, 
no restriction on the masses, concentrations, and force constants will be considered. Consequently, the transport 
coefficients are nonlinear functions of both the shear rate and the thermal gradient and of the parameters of the 
mixture. Obviously, our results reduce to those obtained from the BGK model in the case of particles mechanically 
equivalent Further, when one considers a binary mixture and takes the tracer limit (molar fraction of one of the 
components much smaller than 1), previous results derived from the GK model are also recovered ||. 

The plan of the paper is as follows. Section || starts with a brief survey of the GK model for an A^-component 
mixture. Then, the nolinear Couette flow is introduced. In the same way as in the single gas case, we propose a 
solution characterized by constant pressure and linear velocity and para boli c temperature profiles with respect to a 



scaled variable. The consistency of the assumed profiles is fulfilled in Sec. |II| from a formal solution of the GK model. 
The consistency condition for the temperature yields a system of N implicit equations which give the shear rate 
dependence of the partial temperatures Tj of each species. For nonzero shear rates, the Tj's are different, so that the 
kinetic energy is not equally distributed among the different species. This effect is generic for multicomponent systems 
in far from equilibrium states and is consistent with other studies of mixtures out of equilibrium ]l0|-|l2]] . Apart from 
the partial temperatures, the pressure tensor and the heat flux vector are also explicitly determined. From these 
fluxes the main transport coefficients of the problem can be identified. Specifically, there are five relevant transport 
coefficients: the generalized shear viscosity function, two viscometric functions, the generalized thermal conductivity 
function, and a coefficient measuring cross effects in the heat flux (heat flux along the x direction due to a thermal 
gradient along the y direction). The expressions of these functions are written in Sec. The particular case of a 
binary mixture (N = 2) is considered in Sec. [v|, where the shear rate dependence of the above trans por t coefficients 
is illustrated for several values of the parameters of the mixture. Finally, we close the paper in Sec. |V| with a brief 
discussion. 



II. THE KINETIC MODEL AND THE PROBLEM 

We consider an iV-component mixture. In the low-density regime, the one-particle velocity distribution function 
/i(r, v; t) of species i [i = 1, . . . , N) obeys the Boltzmann equation || 

^1 • v • V/, = Kulfi, fi] + X K iAfi, M> (1) 

where Kij[fi, fj] is the nonlinear Boltzmann collision operator. The first and second terms on the right-hand side 
represent self- and cross-collisions, respectively. They conserve the number of particles of each species, the total 
momentum and the total energy The local number density ni and mean velocity of species i are defined as 



These quantities define the total number density n — ^ i n,; and the flow velocity 

N 



n t = / rfv/, , (2) 



u; = — / dvvfi . (3) 



I 

u=-^p ( u,, (4) 



»=i 



where pi = miiii and p — ^2 i Pi 1 mi being the mass of a particle of species i. It is also convenient to define a 
temperature Ti for each species through 
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^nik B Ti = J dv(v - u^) 2 f t 



(5) 



which is a measure of its mean kinetic energy per particle. Here, k B is the Boltzmann constant. From these partial 
temperatures, the temperature T of the mixture (which is the relevant one at a hydrodynamic level) is given by 



A'* 



ik B T = } 



nik B Ti + ^pi(ui - u) 2 



(6) 



The corresponding balance equations associated with n$ , u and T can be easily obtained from the set of Boltzmann 
equations ([!]): 



D t m + UiW ■ u 



V-Ji 



(7) 



Au + p^VP = 0, 



(8) 



D t T - - V 

n m.j 



(V • q + P : Vu) = 0. 



3nk B T 

In the above equations, D t = 9 t + u • V, and we have introduced the dissipative fluxes of mass 

Ji = mi / dVV/< ' 

momentum (pressure tensor) 



and energy (heat flux) 



N N 

P = ^2 rrii / dvWfi =J^P 

A — 1 <J A — 1 



TV - N 

q = Ey ^ 2v /. = E*. 



(9) 



(10) 



(11) 



(12) 



where V = v — u is the peculiar velocity. These fluxes define the relevant transport coefficients of the mixture. 

In general, the set of N nonlinear Boltzmann equations (]l|) cannot be solved, especially in far from equilibrium 
situations. This is basically due to the complex mathematical structure of the Boltzmann collision operators. In order 
to overcome such difficulties one can resort to kinetic models in which case the Boltzmann collision operators are 
replaced by simpler terms that preserve their main physical properties. Here we consider the well-known Gross-Krook 
(GK) model ||, where the Boltzmann terms Kfj[fi, fj] in Eqs. (|l|) are replaced by simple relaxation terms of the form 



where Vij is an effective collision frequency and the reference distribution functions are given by 



2-KkBTi-. 



.3/2 



exp 



2k B Ti 



Here, we have introduced the fields 



rriiUi + m,-u 



3 U J 



m,i + m.j 



(13) 



(14) 



(15) 



T = T- 4- 

{rrii + m 3 Y 



(16) 
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The above expressions are obtained by requiring that momentum and energy moments of K^ K be the same as 
those of the Boltzmann operator for Maxwell molecules (i.e., particles interacting through a potential of the form 
&ij = Kiji" - ^). This allows one to identify as 



= An j 



m 



in 

miirij 



1/2 



(17) 



where A — Air x 0.422. Although we are considering Maxwell molecules, we can assign an effective diameter cr^ to the 
interaction between particles of species i and j. Dimensional analysis allows one to interpret the ratio («y / Hjj) 1 ^ 4 as 
the size ratio &ij /<Jjj- On the other hand, it must be remarked that the results derived in this paper could in principle 
be extended to more general potentials. 

We now describe the problem we are interested in. Let us assume that the iV-component mixture is enclosed 
between two parallel plates in relative motion and kept at different temperatures. Let the cc-axis be parallel to the 
direction of motion and the y-axis be normal to the plates. We want to analyze a steady state with velocity (u x ) 
and temperature (T) gradients along the y direction (steady Couette flow). These boundary conditions give rise to 
heat and momentum transport, the relevant fluxes being the pressure tensor and the heat flux. Our goal is to get 
these transport properties in terms of the imposed gradients and the parameters of the mixture, namely, masses, 
concentrations, and "sizes" of the particles. In the geometry of the steady Couette flow and in the context of the GK 
model, the set of N coupled equations for the mixture reads 

d 



-v. 



N 

/. ■ $>«/«, (is) 



where 



JV 

Vj = ^ V 4 j = Vu + ^ V ij ( 19 ) 

is the total collision frequency for a particle of species i. In the same way as in previous works, here we are interested 
in studying the properties of the mixture in the bulk region far away from the plates. Thus, instead of introducing 
the appropriate boundary conditions for Eqs. (|lj), we guess the spatial dependence of the hydrodynamic fields and 
then we verify their consistency. In this sense, we expect to describe transport phenomena in the bulk of the system 
by looking for a consistent solution regardless of the actual features of the boundaries. 

In the case of all the species being mechanically identical, = u, Tj = T and the GK model reduces to the 
well-known (closed) BGK equation || for the total distribution function / = ^\ /»■ I 11 this case, the BGK model 
admits a consistent solution characterized by a constant pressure, and linear velocity u x and parabolic temperature T 
profiles with respect to a scaled variable . This variable takes into account the fact that the rate at which collisions 
take place is nonuniform. The above solution is not restricted to small velocity and temperature gradients and thus 
provides a rare detailed description of nonlinear heat and momentum transport far from equilibrium. Guided by these 
results for the one-component case, a similar solution is proposed here for the mixture. First, on physical grounds one 
expects that the steady Couette flow for the mixture is characterized by (a) constant molar fractions Xi = ni/n, (b) 
the absence of mutual diffusion, i.e., U; = u, and (c) constant temperature ratios Xi = Ti/T. In addition, we assume 
that the partial pressure pi — riiksTi, the velocity Ui_ x , and the partial temperature Tj have the forms: 

Pi = riiksTi — const , (20) 
iH x = <H = const , (21) 



1 d 



My) dy 



Ti = -^7i(<»i) = const • ( 22 ) 



The dimensionless functions 7^ must be explicitly determined by consistency. The conditions = u and Xi = const 
imply, respectively, that 
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v i a i = v j a : j, (23) 

™*^f = , 24) 

X* Xj 

where use has been made of the fact that the ratios Vi/vj are constant. According to Eqs. ( p3| ) and fl24|), only one of 
the shear rates en and only one of the parameters 74 (i = 1, . . . , AT) is independent. The temperature ratios \i depend 
on the shear rate and on the parameters of the mixture. Obviously, %i = 1 if (i) the mixture is at equilibrium or (ii) 
the particles are mechanically equivalent (to.; = to and mj = k). The explicit calculation of the temperature ratios is 
one of the main objectives of this paper. 

To verify the consistency of the assumed solution it is necessary to prove that 



f dv{l,V,V a }/ i = {n i ,0,3^ 

J L m i 



(25) 



The fulfillment of these conditions will be carried out in the next Section. 



III. CONSISTENCY OF THE SOLUTION. CALCULATION OF THE FLUXES 



In this Section we show that the relations (|25|) are verified and obtain the pressure tensor and the heat fux. To do 
that, let us consider the formal solution to Eq. ( |l8| ) given by 



Notice that when the operator [1 + (v y / i>i){d / dy)]^ 1 acts on the hydrodynamic quantities the expressions (20)-d2S 
are assumed to hold. By using this formal solution in the Appendix we show that the conditions for the density and 
the flow velocity are satisfied with independence of the value of the parameters Tj in Eq. ( |22] ) for the temperature 
Tj (y) . A similar analysis for the temperature condition shows that a solution exists only for a particular choice of 7^ 
given implicitly by the equation: 

3 N 

2 E t 1 - *i) = E ^Xti K^OtoO - 7« [2F 2 ( 7 ii) + 3Fi( 7y )]} . (27) 



Here, we have introduced the quantities jij — Xijlii 



miVcin 



'(mj + TOj) 2 Vxi 

In addition, the functions F r (jij) are defined as 



Mnj) = [ 4-ia ) Mm) (29) 



with 

9 /*°° 

Fo(7ii) = — / rf«e- t2 / 2 if (2 7 ^ 1/ V/ 2 ), (30) 

i^o being the zeroth-order modified Bessel function Jll|. The set of TV coupled equations (^7j) must be numerically 
solved subject to the constraints (Q and 

N 

^TxiXi = l. (31) 

i=l 
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From a mathematical point of view, the problem is now well posed since we have N independent equations (E7|) for 
a set of N unknowns, say for instance {71, xi, X2, ' " ■> Xn-i}- The solution of these N implicit equations gives the 
above quantities as functions of the shear rate and the parameters of the mixture (masses, concentrations, and sizes). 
This completes the confirmation of the consistency conditions for the hydrodynamic fields. 

Once the functions 7 $ and Xi are known, one can derive explicit expressions for the fluxes. In the case of the Couette 



flow, the mass flux Jj vanishes, and so the relevant fluxes are the pressure tensor (|llj) and the heat flux ( |12|) . We 



begin with the nonzero elements of the pressure tensor. In the Appendix it is shown that they can be written as 

N N 

P xx = nk B TY J ^ 1 Y. "aW i 1 + 4 ^ t Fl (7« ) + F2 (lij )] } , (32) 
i=i '« j=i 

N N 

P yy = nW^^ J2 "Mi i 1 - [Ktnu) + 2f a(7«)]} » (33) 
i=i j=i 

N N 

P zz = nk B T Y^^rY. u VX<j [1 - 2 7y *l {lH )] , (34) 

i=l Vt j=l 

N x N du 
p xy = -nk B T Y ~T" Y "vXijFo {lij ) ■ (35) 
i=i * j=i y 

Notice that the pressure tensor does not depend on the thermal gradient. With respect to the heat flux, its nonzero 
components are q x and q y . The fact that q x ^ is a direct consequence of the presence of the shear flow. These 
components are also evaluated in the Appendix and the result is 



5 N x 2 N 

% = -^ nk B T Y~^r^ 51 Ui jXij { 3i? l(7y) + 2 ^2(7^ 



hm l vf . . 
i=i 1 j=i 

+2a 2 [4^ 4 ( 7tf ) + 4F 3 ( 7<j ) + F 2 ( 7 „ )] } — , (36) 

N 2 N 

q x = nk 2 B Tj2 a *^2 Y Vi iXij {^(jij) + 2F 3 ( 7 ij) 
i=l * * 4=1 

9T 

+2a 2 [4F 5 ( 7ij ) + 8^ 4 ( 7iJ ) + 5F 3 ( 7<i ) + ( 7tf )] } — . (37) 

From Eqs. ([35j) and d36|) , it is easy to prove the identity d y q y = —P xy d y u x , which is in fact required by the energy 
conservation equation (ft) to support the forms for the hydrodynamic profiles. 

IV. GENERALIZED TRANSPORT COEFFICIENTS 

It is useful to summarize the above results for the fluxes in terms of appropriate scalar transport coefficients. To do 
that, let us reduce the velocity gradient d y u x with respect to an effective collision frequency v. Since our description 
applies for arbitrary mass, concentration and size ratios, we choose for simplicity the average collision frequency 

1 N 

and define the (dimensionless) shear rate 

1 d 

a = ———u T = const (39) 
vdy v ' 



G 



as the relevant nonequilibrium parameter. Obviously, a, = (y/v{)a. 

Momentum transport is typically characterized by three rheological functions. The xy element for shear stressess 
can be written in the form of a generalized Newton's viscosity law 

P*v = -VoF v (a)^, (40) 

where 

JV 



i=i 



r ]0 =nk B Tj2- (41) 



is the Navier-Stokes shear viscosity of the mixture Q and the explicit expression of F v (a) can be obtained from Eq. 



N \ N N 



F i(°H£S E^fE^(7ii). (42) 



Vi / ' — ' v, 

=1 / 1=1 8 j=l 



When a = 0, F v = 1, so that this function describes the nonlinear rheological effects associated with the transport 
of the x component of the momentum along the y direction. Normal stresses are measured by the dimensionless 
viscometric functions ^\^{o)'- 



* x (o) = ±yy ; xx , (43) 

* a ( ) = F "~ p " , (44) 

per 

where p = nksT. The explicit expressions for the viscometric functions can be easily obtained by substituting Eqs. 
©-© into Eqs. @ and @. The result is 

= ~a E — E ^X«7y [3Ji(7tf ) + 4F 2 ( 7ii )] , (45) 

° i=l 3=1 

N JV 

^W^EvE VijXiiUjUHjY (46) 
° i=l J 3=1 

The heat flux introduces two relevant transport coefficients, F\ and $, defined respectively as 

q y = -\ F x (a) — , (47) 
<9y 



<9T 

= -A $(a)a— . (48) 

Here, Ao is the Navier-Stokes thermal conductivity coefficient of the mixture given by j^] 

r JV 

Ao = >bE^- ( 49 ) 
%=i 

Equation ( p7| ) can be interpreted as a generalization of Fourier's law with the thermal conductivity modified by the 
rheological quantity F\(a) 1 a nonlinear function of the shear rate and the parameters of the mixture. Equation ( |4"s| ) 
provides information on the anisotropy induced by the Couette flow in the transport of energy. It gives a heat flux 
along the x axis due to a thermal gradient along the y axis. This effect is absent in the Navier-Stokes regime, being 
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of first order in both the shear rate and the temperature gradient. Additional effects are described by the function 
<I>(a). The explicit expressions of Fx and $ can be derived from Eqs. (|3^) and (^) as 

(N \ _1 N 2 N 

2 2 ^ 

+2^ [4F 4 ( 7 y ) + 4F 3 ( 7o - ) + F 2 ( 7tf )] | , (50) 

*<«) - -\> IE E g| 5>«„ f^h«) + »»<■*> 

\i=l / i— 1 1 j — 1 

2 2 ^ 

+2^- [4F B ( 7li ) + + 5F 3 ( 7ij ) + F a ( 7ii )]| . (51) 

Before closing this Section, let us explicitly write the corresponding expressions of the generalized transport coef- 
ficients in the simple case of particles mechanically equivalent. In this limit, rrii = to, Kij = k, v = An(2K/m) 1 ^ 2 , 
7 i = 7 , and Xi — 1- Consequently, Eqs. (|||), (j45|), (f46|), (^0|), and ( |5l"| ) reduce, respectively, to 

^(o) = F ( 7 ), (52) 

* i(a) = ~-y [3Fi( 7 ) + 4F 2 ( 7 )] , (53) 
or 

* 2 (a) = 4^ F 2(7), (54) 

F x (a) = i {3Fx( 7 ) + 2F 2 ( 7 ) + 2a 2 [4F 4 ( 7 ) + 4F 3 ( 7 ) + F 2 ( 7 )]} , (55) 

*(o) = -| {5F 2 ( 7 ) + 2F 3 ( 7 ) + 2a 2 [4F 5 ( 7 ) + 8F 4 ( 7 ) + 5F 3 ( 7 ) + F 2 ( 7 )]} . (56) 
All these expressions coincide with those previously derived from the BGK equation I^Jl^fl. 

V. AN ILLUSTRATIVE EXAMPLE: A BINARY MIXTURE 

Thus far, all the results hold for an TV-component mixture. For the sake of illustration, in this Section we will 
consider the case of a binary mixture (N = 2). Therefore, the corresponding set of two implicit equations ( p7| ) give 7 i 
and xi as functions of a and the parameters of the mixture, namely, the mass ratio /i = toi/to 2 , the concentration ratio 
S = X1/X2, and the size ratios u>n = Kxx/kxi, and w 22 = k 22 /ki 2 . Although the natural nonequilibrium parameter is 
the shear rate a, from a practical point of view it is convenient to take 7 i as an independent variable and express all 
the unknowns in terms of this parameter. In this sense, xi is defined by the implicit equation 

vi\ 2 ELi viiXuFidii) _ 3^i2(l - X12) +_2_SLi ^uXulu [2f 2 (7i») + 3J\( 7 n)] 
1/2 ) ELi ^X2iFi( 72l ) ~ 3^ 2 i(l - X21) +2Ei=i ^2iX2i72i [2F 2 ( 72 i) + 3Fi(7ai)] ' 



and the shear rate a is given by 



,2 



1/^x2 3i/i 2 (1 -X12) + 2ELi^hXh7h [2F 2 ( 7 ii) + 3F 1 ( 7li )] 



2 v ^ y Ei=i ^iiXiiFi( 7 i 



(58) 



In the above equations, 72 and xi are related to 7 i and xi through Eqs. (|24j ) and (31). In general, Eq. (|57| ) must be 
solved numerically. However, 7l and xi can be explicitly derived in some limiting cases. For instance, in the limit of 
small shear rates 7 i — > -y^a 2 , and %i — ► 1 + x^a 2 , where 
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(0) 

7i = 



v \ vi + 8v 2 
V\ ) \iv\ + 5v 2 



(59) 



Xi 



(0) _ l(/i 2 -l)(l+/i)(l+a) 



Vl2 



fXVl + Sv 2 ' 



(60) 



The temperature ratio T\/T 2 is a measure of the lack of equidistribution of the kinetic energy. This quantity is 
plotted in Fig. [l] as a function of the shear rate for w\\ = w 22 = 1, a concentration ratio 8 = 2 and two different values 
of the mass ratio: /i = 0.5 and /i = 2. We observe that the kinetic energy per particle of the solvent component 1 is 
larger (smaller) than that of the solute component 2 if the former is heavier (ligther) than the latter. The fact that 
Ti T 2 for arbitrary shear rates indicates that the standard one-temperature Chapman-Enskog theory |g] fails, even 
for not very disparate masses. In this case, the so-called two-fluid theory must be used fl5|| . On the other hand, the 
monotonic behavior for T\/T 2 found in the steady Couette flow contrasts with the one observed in the uniform shear 
flow [prf . In the latter problem, for given values of the parameters of the mixture, there exists a particular nonzero 
value of the shear rate a such that T\ = T 2 . 




FIG. 1. Shear-rate dependence of the temperature ratio T1/T2 for wn = W22 = 1,8 = 2, and two values of the mass ratio fi: 
(a) fj, = 2, and (b) [s, = 0.5. 



The viscosity function F v of a binary mixture under Couette flow is defined by Eq. (j42j). For small shear rates this 



function behaves as F„ 



l + F^a 2 , where 



p(0) 
v 









Mx<°> ( 




V2 ) 







X2 
V 2 



(0) 



2^(1 + *) 
v 2 



5_ 

M 



I87I 



(o) 



18 



7i 



(0) 



(61) 



with M = /i/(l + /i) 2 . Depending on the values of the parameters of the mixture, the super-Burnett contribution 

can be positive or negative, so for small rates the viscosity function can increase or decrease as the shear rate 
increases. Nevertheless, outside this small region, F v monotonically decreases as the shear rate increases, regardless 
of the values of the parameters of the mixture as shown in Fig. ^[ 
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2 

a 

FIG. 2. Shear-rate dependence of the viscosity function F v for Wu = W22 = 1,8 = 2, and three values of the mass ratio \i: 
fi = 0.5 (solid line), fi = 1 (dashed line), and [i = 2 (dotted line). 



In addition, at a given value of the shear rate, the viscosity function increases as the mass ratio mx/mz decreases 
when ni/ri2 > 1. The viscometric functions $1,2 are given by Eqs. (45) and (^6|), respectively. In the limit a — > 0, 
one gets 



* 1 (0) = -14(a; 1 +X2/i^ ) 7i°\ 



(62) 



*2(0) = 4 x x +x 2f i-^ 7} 



, ; (0) 



(63) 



Notice that '5 1.2(0) are Burnett transport coefficients. The shear rate dependence of the (dimensionless) viscometric 



functions ^1,2(0) is shown in Figs 
is always negative while the second one is positive 
with the shear rate. On the other hand, for zero shear rates 
inequalities are kept for finite shear rates. 



and ||, respectively, for the same cases as in Fig. |2j. The first viscometric function 
In general, the magnitude of both viscometric functions decreases 
1*1,2(0)1^=0.6 < |*i, 2 (0)U=i < 1*1,2(0)1^=2. These 
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0.0 0.2 0.4 2 0.6 0.8 1.0 

a 

FIG. 3. Shear-rate dependence of the first viscometric function — *i for wu = W22 = 1,8 — 2, and three values of the mass 
ratio fi: n = 0.5 (solid line), fi = 1 (dashed line), and fi = 2 (dotted line). 




a 

FIG. 4. Shear-rate dependence of the second viscometric function ^2 for wu = W22 = 1,8 = 2, and three values of the mass 
ratio fj>: /i = 0.5 (solid line), fi = 1 (dashed line), and fi = 2 (dotted line). 

The transport of energy is characterized by the coefficients Fa and For small shear rates, F\ — > 1 + F^a 2 and 
$^$(0),with 



,(o) 



/ Xi 



X2 



\miVl VT12V2 



x\ 
m\v\ 



+ ■ 



X2 

m 2 v 2 



5} 252 ^^ (0) 36 v 2 
M 



(64) 
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a 



FIG. 5. Shear-rate dependence of the generalized thermal conductivity function F\ for Wu = W22 = 1, S = 2, and three 
values of the mass ratio (i: fj. = 0.5 (solid line), [i — 1 (dashed line), and /j, = 2 (dotted line). 



It is apparent that the coefficient presents a complex dependence on the parameters of the mixture. In general, 
and in accordance with Fig. |5[ for finite values of the shear rate the generalized thermal conductivity coefficient Fx 
decreases as the shear rate increases with independence of the values of the parameters of the mixture. It is interesting 
to remark that the curves of /x = 0.5 and fi — 1 are practically the same in the range of shear rates considered. In Fig. 
^ we plot the shear rate dependence of the reduced coefficient $(a). We observe that <E> is negative and its magnitude 
decreases with the shear rate. Further, in the range of shear rates studied, this coefficient is practically independent of 
the mass ratio. This is rather surprising if one takes into account that this is a coefficient measuring complex coupling 
effects between both gradients and a more significant influence of the parameters of the mixture should be expected. 
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FIG. 6. Shear-rate dependence of the cross coefficient — <3> for wu 
/i = 0.5 (solid line), fi = 1 (dashed line), and fi = 2 (dotted line). 



1,5 = 2, and three values of the mass ratio /j,: 



Finally, let us note that in the tracer limit (n\jni 
previously obtained from the GK model @. 



0), all the results of the binary mixture reduce to those 



VI. DISCUSSION 



In this paper we have considered a gas mixture of N species in a steady nonequilibrium state involving combined 
heat and momentum transport. This state is usually referred to as the steady planar Couette flow. The system is a gas 
mixture enclosed between two relatively moving parallel plates which are kept in general at different temperatures. The 
shear rate and the thermal gradient are the two nonequilibrium parameters of the problem. Under these conditions, 
the density, flow velocity, and the temperature are nonhomogeneous, thus making the analysis of this problem more 
complicated than in other nonequilibrium situations. Further, the mixture is not restricted to specific values of the 
mass ratios, concentration ratios and size ratios. As a consequence, the transport properties of the system are nonlinear 
functions of the strengths of the gradients and of the parameters of the mixture. In order to offer a description as 
detailed as possible, we have used a kinetic model of the Boltzmann equation. Specifically, we have considered the 
GK model for mixtures, whose reliability has been proved in the past in several nonequilibrium problems JT^Jl^ ]. The 
explicit determination of the transport properties of the gas from an exact solution of the GK model has been the 
goal of this paper. It is important to remark that progress has been possible here due to the previous results derived 
in the limit cases of mechanically equivalent particles H and tracer particles || . 

Guided by the results obtained for a single gas M, we have found a consistent solution of the GK equation 
characterized by a uniform presure, and linear velocity and parabolic temperature profiles with respect to a certain 
scaled variable (distance times the local collision frequency). The consistency condition for the partial temperatures 
Tj leads to a set of N coupled transcendent equations, whose (numerical) solution gives Ti in terms of the shear rate 
and the parameters of the mixture. Once the TVs are known, we explicitly evaluate the pressure tensor and the heat 
flux. The nonzero elements of the pressure tensor define three relevant transport coefficients: the generalized shear 



viscosity F n , Eqs. (M) and (|42|), the first viscometric function Eqs. (J43|) and (|45|), and the second viscometric 
function Vt^, Eqs. ((44|) and (|46|). While the viscosity function measures shear stresses, the two viscometric functions 
provide information on normal stresses. The heat flux turns out to be exactly proportional to the thermal gradient, 
with independence of its magnitude, but it defines two shear rate dependent transport coefficients: the generalized 
thermal conductivity Fx, Eqs. ( [i"7|) and (|5(]), and a generalized cross coefficient Eqs. (fj§| ) and (fil|), which accounts 
for the transport of energy along the direction normal to the thermal gradient. For the sake of illustration, a binary 
mixture has been considered. The results show that, except perhaps in a narrow region small of shear rates, the 
magnitude of these coefficients decreases for finite shear fields. Consequently, the main effect of the shear flow is to 
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inhibit the transport of momentum and energy across the system. In general, this inhibition is more significant when 
the excess component is heavier than the defect component in the case of the coefficients and F\ , while the opposite 
happens in the case of the two viscometric functions. Concerning the cross coefficient 'f (which is negative), the results 
indicate that in the range of shear rates considered is practically independent of the values of the parameters of the 
mixture. 

The problem studied here is worthwhile by itself. Furthermore, it can be taken as a starting point to analyze 
mutual diffusion under steady Couette flow. In the same way as done in the tracer limit |l^,|l^] , the idea is to carry 
out a perturbation expansion in powers of the molar fraction gradient around the solution found in this paper. In 
the first order of the expansion, a generalized Fick's law is expected, where a shear rate dependent mutual diffusion 
tensor can be identified. On the other hand, it is apparent that the results reported here can also be of relevance 
in connection with computer simulations. In the context of molecular dynamics simulations, it is difficult to achieve 
large shear rates in the bulk region to clearly observe nonlinear effects in the transport of momentum and energy. One 
possibility to overcome such problems inherent to molecular dynamics in the low-density regime is to use the direct 
simulation Monte Carlo method [^C|, which has been shown to be fruitful in several problems. We hope that the 
results derived here for the nonlinear transport coefficients may stimulate the performance of computer simulations 
to check the reliability of our predictions. Recent comparison made for a single gas under Couette flow supports the 
accuracy of the BGK results. 
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APPENDIX A: VELOCITY MOMENTS OF THE DISTRIBUTION FUNCTIONS 

In this Appendix we evaluate the velocity integrals for the consistency conditions and for the momentum and heat 
fluxes. They correspond to the low order moments of the distribution functions fi. Let us consider the formal solution 
to Eq. (|lj) given by 



N 
3 = 1 



(Al) 



where 



A,„ = (l + ^ 

Vi ay 



fe=0 



(A2) 



We define the velocity integrals 



dvV^V y k2 V z k3 A i:j . 



(A3) 



From a formal point of view, the moments defined in Eq. (A3) can be directly obtained from the results derived in 
the one-component case when one exploits the equiva lenc e between the hydrodynamic profiles (20) (|2^) with 

those given in the single gas. Consequently, the moments (A3) can be determined from comparison with those of 
the distribution function / of the one-component gas by making the changes v —> Vij, v~ x d y u x — > <Zj, T — > T^ , and 
{v~ l d y ) 2 T = — (2m//cs)7 — > — (2m,i/kB)^ij- Taking into account this equivalence, we can explicitly write the first few 
moments of Ay . Here, we display the moments with degree k\ + ki + k% < 3 appearing in the calculations performed 
along the main text. They are given by [Q,[l4| 



14 



^0,0,0 ~ n ii 



(A4) 



(A5) 



(A6) 



'0,2,0 



(A7) 



(A8) 



^"1,1,0 ^ -^Ofli; 



TO, 



(A9) 



m 0,l,2 



9 * 1 r, 



(A10) 



(All) 



•m~Vi ay 



1,2,0 



(A12) 



M« )2 = 2^piF,,-T, 



mfui dy 



(A13) 



M ^ = _^|_J1 [ Fl + 2a ? (4F 4 + 4F 3 + F a )] , 
' ' mfvi ay 



(A14) 



M ^ n = 2 n 1 ksl\ 1 ^ + ^ + g ^ + ^ + p ^ a ._H_ T .. 



'3,0,0 — mfvi 



(A15) 



In these expressions, 7^ = Xijlu Xij 1S defined in Eq. (28) and 

d 



(A16) 



with F (7i 3 ') given by Eq. @. 

Th e co nsistency conditions for the density and the flow velocity are easily seen to be verified by using Eqs. (A4) 
and (A5), respectively. The condition for the temperature on the other hand requires that 



N 

E?( M 2'o,0 + ^2,0 + ^fo, 2 )=3 



nik B Ti 



(A17) 



Substitution of Eqs. (|A6|) (JAq) into Eqs. (A17) yields the set of coupled equations (|27j). The momentum and heat 
fluxes can be evaluated from the remaining velocity moments. The pressure tensor can be written as 
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N . N 



(A18) 



where the nonzero elements of Pj are given by 



JV 

{Pi,xx'i Pi,yy\ Pi,zz] Pi,xy} = ^ m i j ^2^0,0' -^0^2,0! -^ofo,2i ^1,1,0 \ 



The heat flux can be obtained in a similar way. It can be also written as 

N „ AT 



(A19) 



(A20) 



where 



Qy,i 



N 



(A21) 



AT 



(A22) 



The explicit expressions for the fluxes and transport coefficients appearing in the text can be easily derived from 
Eqs. (|A1^)-(|A22|). 
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